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Abstract In electromagnetism a current along a wire tightly wound on a torus makes 
a solenoid whose magnetic field is confined within the torus. In Einstein's gravity we give 
a corresponding solution in which a current of matter moves up on the inside of a toroidal 
shell and down on the outside, rolling around the torus by the short way. The metric is static 
outside the torus but stationary inside with the gravomagnetic field confined inside the torus, 
running around it by the long way. This exact solution of Einstein's equations is found by 
fitting Bonnor's solution for the metric of a light beam, which gives the required toroidal 
gravomagnetic field inside the torus, to the general Weyl static external metric in toroidal 
coordinates, which we develop. We deduce the matter tensor on the torus and find when it 
obeys the energy conditions. 

We also give the equipotential shells that generate the simple Bach- Weyl metric externally 

and find which shells obey the energy conditions. 

1 Introduction 

We study this problem firstly to illustrate the power of thinking about the gravomagnetic 

field as analogous to the magnetic field even when gravity is strong, and secondly to show that 

intuition gained from studies of cylindrically symmetric space-times can often be justified. 

Cylindrical systems have infinite total mass and their metrics are not flat at infinity. Often 

they are not even limits of finite equilibria as some parameter tends to infinity. However, 

when a cylinder of finite length is bent around into circle with the ends joined to make a 

torus, it gives a system that can sometimes be solved even for strong fields. When the original 

cylinder was rotating about its axis the corresponding torus will be rolling about its central 

circle. Here we study such systems when all the mass resides on a shell which is either a 

torus or in the static Bach- Weyl metric an equipotential. 
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Toroidal solutions of Einstein's equations have been considered before. As early as 1922 
Bach and Weyl [3] gave the simplest static solution which has been further elucidated by 
Hoenselaers [7J, and studied by Semerak et al |10j . The general exterior solution for the 
static metric in toroidal coordinates was given by Frolov et al [5] in an investigation of 
cosmic strings but they did not perform the integration that is necessary to get all the metric 
coefficients explicitly. Here we shall need that full solution for our external field. We also give 
explicit examples of equipotential shells that generate the static Bach- Weyl metric externally 
and find out when such shells fail to obey the energy conditions. However this paper is 
primarily devoted to the rolling tori which have toroidal gravomagnetic fields inside a rolling 
matter shell but are externally static. 

The 1966 edition of the Classical theory of Fields by Landau and Lifshitz [H] gives Ein- 
stein's equations for general stationery metrics in a form that has strong analogies with 
Maxwell's electrodynamics. The technique identifies the points of space that lie along the 
time-like Killing vector so it does not extend continuously inside ergospheres where the Killing 
vector becomes space-like. Later Geroch [B] and others exploited the special properties of 
these equations in developments that led to the generating techniques for new solutions. We 
write the metric in the form 

ds 2 = f(dt + A k dx k ) 2 - j M dx h dx l = g^dx^dx" , (1.1) 

where k and I run from 1 to 3. Since the metric is stationary £, A k and -fki are all independent 
of t but in general they depend on the x k . We work in the positive definite three dimensional 
metric of space, 7/y . It is not a cross-section of the four metric by any surface, nevertheless we 
may define its Christoffel symbols and the corresponding three-dimensional Ricci tensor 
of this gamma space, use commas to denote ordinary derivatives and semicolons to 

denote covariant derivatives in gamma space. The Ricci tensor of space-time will be denoted 
by Rjiv\ its indices are raised and lowered by g^ while the indices on Ak and P kl are raised 
and lowered by the gamma metric. One may show that g kl = —j kl ,k,l = 1,2,3 and that 
the determinants of the metrics are related by —g = ^7: In gamma-space we define the 
alternating tensor rf^ k = e^ k /^/j where epsilon is the alternating symbol which is unity 
when i,j,k are an even permutation of 1,2,3, minus one for an odd permutation and zero 
otherwise, 77^ = i/je^k. The divergence and curl are defined in gamma-space by 

dWE = (l/^)d k (^E k ) ; (curli^ = {l/^)e^ k djE k ; (grad^) fc = d k <f>; (1.2) 
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so curl grad and div curl are both zero. We define the gravomagnetic induction 3 by 



B = curl .4., (1.3) 
where Ak is the vector potential defined in the metric (1.1). Clearly diviE? = so 3 carries 



the gravomagnetic flux. Landau and Lifshitz rewrite Einstein's equations in gamma space; 
rewriting their equations in our notation we have with k = 8ttG/c , 

£divgrad£ + i £ 4 2? 2 = R 00 = K (T 00 - 1 <? 00 T). (1.4) 

Henceforth we use units with c = 1 and G = 1. If we now define a field intensity vector 
T-L = then their second equation reads 

(curm) fc = -2<T fc = -2 K J k . (1.5) 

Notice a strong resemblance of this strong field equation to Maxwell's electrodynamic equa- 
tion curlH = 4-7rj. In both cases the current has no divergence however in general 'H. has 
a divergence while IS does not. Clearly H is the gradient of a scalar whenever J is zero. 
The T-L k are the spatial components of the twist vector T] tJ,VUT ^ v D T ^ CT where D is the covariant 
derivative in the space-time g^ v . The last Einstein equation is 

P kl + i t 2 {~i kl B 2 - B k B l ) - C 1 ^' 1 = R kl = <T kl - i g kl T). (1.6) 



Einstein's equations (1.4), (1.5) and (1.6) do not mention .4. itself but only 3 = curl.4., 
so their solution for A. is arbitrary up to a spatial gradient. This gauge transformation 
corresponds to the transformation A.' = A. + grad to (x k ). This corresponds to shifting the 
zero of time by a position-dependent function to. Under such as shift the metric remains of 
the same form with all the metric components independent of the new time. For our toroidal 
problem we want a gravomagnetic field inside the torus to be in the toroidal direction along 
the unit vector n^. As the current is confined to the toroidal shell, Rq = inside, so 
curl'H = 0. Hence 

= U = -8/grad^, (1.7) 

where I is a constant which represents the total matter flux moving around of the torus by the 
short way. In the electromagnetic analogue the —8 would be replaced by a 2 but there we have 
Vx H = 4-7rj whereas Einstein's equation has — 16-7T in place of 4tt. We might have proceeded 



by inserting expression (1.7) into equation (1.5) and then attempt to solve it for £ within 
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the torus where Rqq = 0. However empty space solutions with solely toroidal gravomagnetic 
fields are already available. One is obtained by boosting the Levi-Civita solution along the its 
axis by a Lorentz transformation which gives a current along that axis and a gravomagnetic 
field around it. However that solution is unnecessarily cumbersome compared with Bonnor's 
jl] beautiful exterior metric for the gravitational field of a cylindrical light beam. In this ds 2 
takes the form 

ds 2 = dt 2 -(dx 2 +dy 2 +dz 2 )+(F-l){dt-dz) 2 = F[dt-{l-F- 1 )dz] 2 -{dR 2 +R 2 dip 2 +F- 1 dz 2 ), 

(1.8) 

where outside the beam F is harmonic in the two dimensional x, y space and 

R 2 = x 2 + y 2 ; F = 8/ln (R/a) + C, (1.9) 
and C is a constant. For this metric the gravomagnetic intensity is 

H k = (^curlA) k = {0,F 2 Rd A [(F- l)/F],o} = {0,81,0}, (1.10) 



so, cf.(1.7), I or if we like Ic, is the matter flux. As Bonnor shows his source for this external 
solution is a cylinder of null dust travelling in the z direction. It is this matter current that 
generates the toroidal gravomagnetic field 1-L externally. It is this toroidal gravomagnetic 
field that we need INSIDE of our torus. To avoid the confusion generated by using the 
external part of Bonnor's metric for the inside of our torus we shall in what follows use the 
term Bonnor's metric to refer to his external metric. His internal metric plays no part in our 
calculations. Our procedure for finding the metric of a rolling torus is to adopt a version of 



Bonnor's metric inside our torus. We integrate equations (1.4), (1.5) and (1.6) with delta- 



function contributions to Rqq,Rq and R M on the torus, to obtain the junction conditions by 
which the external solution must be fitted to the internal one. We notice that A itself as 
opposed to B(= curU) is not involved in these junction conditions since none of equations 



(1.4), (1.5) and (1.6) involve A itself. In fact the first two junction conditions are closely 
related to the electrodynamic ones that n x E and n B must be continuous where n is 
the unit outward normal and the discontinuities m n ■ E and n x B give surface charges 
and surface currents. In the gravitational case they give a surface mass density and mass 



currents. The junction condition for integrating equation (1.6) is a three-dimensional version 



of Israel's condition relating the discontinuity in the external curvatures of a three-surface 



to the surface stresses and currents. Since from equation (1.10) the internal metric has no 
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gravomagnetic field penetrating the surface of the torus, discontinuities in n x 3 can all 
be catered for by matter currents in the surface just as they are in electromagnetism. To 
clarify this analogy the next section gives the flat-space solution of Maxwell's equations for 
a toroidal solenoid. This also serves to introduce toroidal coordinates that we use for both 
the electrodynamic problem and the Weyl solutions. In our gravitational problem there is 
no gravomagnetic field outside the torus so the metric outside is static. Thus we can use 
the general solution for the static external metric in toroidal coordinates which we develop 
in section 3 and illustrate with a static shell source for the Bach- Weyl metric in section 4. 
The main problem is then reduced that of fitting our general external solution to Bonnor's 
solution which we use for the inside of our torus. This we do in section 5. In section 6 we 
inspect the energy conditions on the surface stresses and surface currents to ensure that they 
are obeyed. 

2 The Toroidal Solenoid in Electrodynamics 

In this section we use the usual Cartesian conventions for our flat-space 3-vectors, rather 
than using covariant and contravariant components. 

2.1 Toroidal coordinates 

While toroidal coordinates are well known, [9] they are rarely used so we have felt it necessary 
to help the readers unfamiliar with them by a brief introduction. Consider two fixed points 
in a plane. The locus of a point having a given ratio of distances r\jri > 1 to those fixed 
points is a circle. By changing the ratio we get a set of coaxial circles. Now consider the 
coordinate system generated by rotating those coaxial circles about the axis that bisects the 
line between the two fixed points. This generates a set of coaxial tori. One coordinate is 
£ = ln(ri/r2) > 0; the other axially symmetric one is — ir < rj < ir, the angle between the radii 
r\ and r?,- The third is the angle < <p < 2n about the axis of symmetry. The cylindrical 
coordinates R, z can be expressed in terms of £ and rj as follows: 



R = h sinh £ 



z = hs'mr] with h 



a 




cosh C — cos rj 



The flat space metric is then 



dR 2 + dz 2 + R 2 df 2 = h 2 (d( 2 + dr, 2 + sinh 2 (dtp 2 ) , 
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Thus h and R are the scale factors in toroidal coordinates, a is the radius of the central line 
torus. On each given torus £ is constant. On the axis R = 0, £ = and £ is also zero at 
infinity. On z = 0, r] = tt when R < a and rj = when R > a. 
We set 

5 = •y/ cosh £ — cos rj, (2-3) 

then for an axially symmetric ip, V 2 ^ becomes in toroidal coordinates 

S 6 



v 2 v> 



a 2 sinh £ 



To separate V 2 tp = we first set ip = SU, then 

! sinhr _ rr „ ! sinr/ x sinhC _ i sinr? 



so after some work, we find that Laplace's equation reads 



S 



s 



V 2 iP = -2 [coth(d ( U + d 2 U + d 2 U + \U] = 0. 



It is useful to denote 



u = cosh ( and v = cos 77, 



then equation (2.6) also reads 



(u - vfl 2 



= V " I' ' {du [{u 2 - l)d u U] + d 2 U + \u) = 0. 



Laplace's equation separates if we write 



U = P(u)E(rj) 



to give 



±{[(u 2 -l)P>]' + lp} = -±d 2 E = l 2 , where P> = g = d u P. 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



(2. 



(2.9) 



(2.10) 



We get thus Legendre's equation for P: 

[{u 2 -l)P'}' - L(L + 1)P = with L = l-\ and Eote ilv . (2.11) 

We note the special solution E = E\j] when I = as well as E =const. We shall be 
interested in u = cosh ( > 1 and L less than an integer by | , so we are interested in Legendre 
functions of u which are not Legendre polynomials. 



The general solution of the P equation in (2.11) is of the form 



P = a l P L (u)+b l Q L (u), 



(2.12) 



but we require solutions that are finite at u = 1 where £ = i.e. on the axis of symmetry 
(and at infinity). Our potential should be symmetrical about z = so only cos(/ry) solutions 
are acceptable. Hence combining the different separable solutions we find that 



z=o 

2.2 The Toroidal solenoid in flat space 

Here we are concerned with Maxwell's electromagnetism in flat space, not with gravomag- 
netism which is the subject of Sections 5 and 6. The current along the wire on the torus 
causes a magnetic field inside. Outside the torus there is no magnetic field so Vx A is zero 
there. 



Figure 1: A toroidal solenoid carrying an electric current where the windings are few and have been 
separated for clarity. In the mathematical problem the current is solely along Vr;. If the solenoid had 
N turns the total current around the pictured torus by the short way would be NI. This quantity is 
called I in the text where the wire is replaced by a continuum. 

However there is a magnetic field passing along the torus so for any loop linking the torus 



CO 



ip = y/u — v U, where 




(2.13) 
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by the short way 

j> A-dl = J B-dS = $ T , (2.14) 

where $t is the total magnetic flux through the torus. Inside the torus the magnetic field is 
in the ip direction. In toroidal coordinates the local unit vectors are 

= KV- ( , n v = /iV- 77 , n v = i?V- <p. 



With both V- B = and VxB, we may write 

^ 21 

= # T V 

/ is a constant current. To find a vector potential A which gives this B we write 

A = AVr], 



(2.15) 



(2.16) 



(2.17) 



V x A will be in the <p direction if A is a function of (, 77 or u, v only. The condition that 

V x A = B is easily found to be 

f)A <?Th 

(2.18) 



d( sinh C ' 

A is thus an integral over £. The arbitrary function of 77 is determined since A must be zero 
on the line torus where £ — > 00. The integral is readily evaluated by writing x = to give 



A 



21a 



In 



x 2 -l 



x 2 — 2ux + 1 V x + 1 



x - 1 



^ln{ " nhC [tanh(l C )r4 
sin r\ y cosh Q — cos rj 1 J J 



(2.19) 



Despite appearances this expression is regular as 77 — > and tt, which is not readily seen from 
the expression in terms of x. 

Outside the solenoid A has to be a gradient since the magnetic field is zero. This is readily 
accomplished by taking A to be the function of 77 that is achieved on the solenoid itself where 
C = (s or x = x s = e^ s . Thus outside 



2Ia 1 
A = — 5— In 

sin 77 



x 2 — 2cosr/x s + 1 \x s + 1 



1 



cos rj- 



For 77 



vr, 



A^2Ia 



In 



1 



(2.20) 



(2.21) 



x s + l (x s + l) 2 _' 

The surface current in the n rj direction around the torus is given in terms of the discontinuity 
of B so 

(2.22) 



4 " J =-^ = -sinhc: Vr? - 



This completes the solution for the magnetic field of a toroidal solenoid. It is a close analogue 
of the gravomagnetic field discussed in Section 5 and 6. 
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3 The general static Weyl metric in toroidal coordinates 



Weyl takes the metric in the form 



e -2^ df 2 _ ^ ^2k ( dz 2 + dR 2j + R 2 d(p 2^ = _ ^ ^2 ^2 + ^2) + R 2^2 

(3.1) 

Then, in empty axially symmetric spaces Einstein's equations give V 2 ip = where V 2 is the 
fiat space operator. Also setting 

D = d R -id z , (3.2) 

we have the Weyl equations 



Dk = \Re^De- 2 ^De- 2 ^. So DkD lnR = (Dip) 2 . 



(3.3) 



From (3.1) we see that £ + ir\ may be obtained by conformal transformation of z + iR so if 
we set 

D* = d c + id v , (3.4) 



equation (3.3) implies via conformal transformation 



D*kD* mi? = (D^f = (d c V>) 2 - (d v *l>) 2 + 218^8^. 



(3.5) 



Now 



InR = [In sinh £ — ln(cosh£ — cost/)] = coth£ — S (sinh £ + i sin if) . (3.6) 



Hence our equations for k become 



sinh 



sin r/ 



d c £^cothC + d v k ^T = (W - (d v i>y 



S 2 



(3.7) 



For regularity we need k = on axis where £ = 0. So & may be found by integrating d^k 
from to C at constant r\. Therefore we eliminate d v k and obtain after a simple calculation 



in terms of u, v defined in (2.7) 
u 2 -l 



d u k = [u 



S 2 



[(u 2 -i)(d u ^) 2 -(i-v 2 )(d v ip) 2 ]+-^(u 2 -i)(i-v 2 )d u ifjd v ifj . (3 
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However our potentials are all of the form ip = SU so putting this form into (3.8) and 
simplifying, 

d u k = -|(1 + uv)U 2 - [v{u 2 - l)d u U - u(l - v 2 )d v U] U 

+ (1 - uv) [(v 2 - l){d u U) 2 - (1 - v 2 )(d v U) 2 ] + 2{u 2 - 1)(1 - v 2 )d u Ud v U. (3.9) 

The general solution for U involves the sum given in (2.13). The general solution for k must 
be obtained by integrating (3.9) which is quadratic in U so that leads to double sums. We 
write 

L = I — g , M = m — 7; and c(n) = cos(n?7). (3.10) 

and after performing the v -differentiations the integral for k is given in terms of Pl and Pm-, 
Legendre functions of u: 

oo oo „ u 

k = J2J2 aia ™ / [P'lP'm{u 2 -1){1-uv)-\P l P m {1+uv)-\ (P L P M )'(u 2 -l)v] du c(l)c{m 

1=0 m=0 ^ 1 

oo oo „ u 

+ Y1 S aiam ta p L p Mm(u 2 -l)+\P L P M (l+m)u\ du [c{l+m+l)-c{l+m -1)] 

1=0 m=0 J 1 

oo oo ~ u 

+ aiClm E" 2 ^L i3 Af" i (n 2 -l) + |PLP A f(/-H«] d ^ [c(Z-m+l)-c(l-m-l)] 

/=0 m=0 ^ 1 

oo oo ^.^ 

+ ^^^^a;a m / \PlPm lm(l — uv) du [c(l + m) — c(l — m)]. (3.11 

/ — n ™ — n J 1 



/=0 m=0 



In Appendix A we show how the indefinite integrals in (3.11) can be evaluated in terms of 



Legendre functions. This results in k being given in the form 

oo oo 

k = I X] X] aiCLm C ( l + m + l ) k lm + c ( l + m ) k l,m +C(l + m- l)k^ 
1=0 m=0 

oo oo 

+ \ a;a m c(Z — m + l)k}_ m + c(/ — "m.)^ _ m + c(Z — m — l)kf^ 



(3.12) 



Z=0 m=0 

where k? are known terms of Legendre functions themselves and their derivatives. Thus tp 



of the form given in (2.13) and k of the form given in (3.12) and the kf given in Appendix 



A constitute the general equatorially-symmetric Weyl solution of Einstein's equations in 
toroidal coordinates . All finite sum solutions are singular on the line toroid at R = a, z = 0. 
We use this general solution in sections 5 and 6 where we discuss the gravitational solenoid 
but before tackling that more complex problem we consider the simplest static toroidal shell 
source. This sheds light on the conicity which is a common feature of both problems but is 
most clearly illustrated without the complications of moving sources. 
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4 The Bach-Weyl metric generated by a static shell 
4.1 The model 



The simplest of the toroidal solutions (2.13) is given by keeping only ao non-zero. Then 



U = a P_ l/2 (u) = a P, 



(4.1) 



which is independent of r\. Henceforward we shall drop the subscript —1/2 and merely write 
P{u) on the understanding that this P is the Legendre function of order —1/2. Substituting 



expression (4.1) into (3.9) all the ^-derivatives vanish and we obtain an expression for dk/du, 
which integrates by parts. Using the recurrence relation for Legendre functions the final 
integral can be put in terms of P and P\ji = P-3/2- The solution for k with the boundary 
condition that k = on axis is then 



k = \al[2P{P,_ -uP)+v [p 2 - (Pi ) 2 ] } 



(4.2) 



Expression (4.5) for k coupled with 



ip = aoy/u — vP, 



(4.3) 



constitute the complete Bach-Weyl metric. If we consider these as solutions everywhere then 
they are singular on the line toroid that forms the circle R = a, z = where u is infinite. 
Indeed using (Jj formula 14.8.14 we find 

o V2 



P{u) 



ITU 



1/2 



ln(8u) 



Thus, for large u 



and at large fixed u, 



7T 



(R - a) 2 + z 1 



u- 



(4.4) 



(4.5) 



(4.6) 



so u is then constant on a torus of small radius a/u = s about the singular circle. The 



behaviour of ip, see (4.5), then shows us, by comparison with both the classical result and 



that of the cylindrical line source in relativity, that the mass per unit length is ^/2oq = 2irfi. 
The metric will therefore suffer from the problems of line sources when the singular line is 
approached too closely. However, when the Bach-Weyl solution is generated as the external 
metric of a massive equipotential shell that obeys the energy conditions, that external part 
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of the metric can not be an unphysical one. The asymptotic form of ip at large r gives the 
total mass m of the system, ij) — > v/ ^°° a so the mass is rh = ™ = ^Jf = \/2ao = 27r/x. While 
this sounds just what we might expect, the proper length of the singular line is e 2 ^2ira — > oo 
since U — > ao an d u 1 ^ 2 ~~ 00 • We shall take a massive shell that lies on an equipotential 
surface ip = V = const of the Bach-Weyl solution. This surface is a toroid in that it has the 
topology of a torus but, as seen in Figure 2, it lacks the circular small cross-section of a true 
torus of constant u. 












-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 

R 



Figure 2: Equipotentials in R, z coordinates for ?/> = 1 the largest toroid crossing the origin, through 
■ijj = 1.01, 1.05, 1.2 and ip — 1-5 the smallest one almost circular. We choose one of these surfaces and 
make it a massive shell with constant internal potential. There is a limit to the mass that can be put 
on it given in Figure 5. 

Inside our shell there is no matter, V 2 ^ = and as tp is constant on the boundary it is 
V everywhere. The equation for k then shows it must be a constant k everywhere inside the 
toroid and space-time is locally flat : 



ds 2 = g^dx»dx v = e- 2V dt 2 - e 2V e 2k (dz 2 + dR 2 ) + R 2 d^ 2 . (4.7) 

Space-time is not globally flat because the axis of symmetry does not touch the locally flat 
space within the toroid, so we cannot use the regularity on the axis to show that k = 0. We 
determine k below. 

The equation of the toroid as seen from outside is, say, z = z(R), defined parametrically 



by (2.1), and, following (4.6), z(R) is defined by the condition that 



V = — = \/h (u — v)P(u) = Const where m = — . (4.8) 
m V * a 



4.2 The junction condition 



The metric on the shell is thus, according to (3.1), with dz = z^dR 



da 2 = h ab d6 a de b = e~ 2V dt 2 -f(l + z 2 R )dR 2 - e 2V R 2 dp 2 , 6° = t,9 1 = R, 9 2 = <p, 
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/ = e v+k and z R = (Orz)^ 



djiip vu 2 -l 
d z tjj sin rj 



-\vP + {I - uv)P' 
\uP+ (u 2 - 1)P' 



(4.9) 



Let z = z{R) be the equation of the toroid as seen from inside. Accordingly, see (4.7), its 
metric is 



da 2 = h ab d9 a d9 b = e~ 2V dt 2 - f(l + z 2 R )dR 2 - e 2V R 2 d V 2 where z R 



■2ud2 „2V d 2t 2 



dz 
dR 



and / = e 



V+k 



Since (4.9) and (4.10) represent the same hypersurface we must have 



e 2fc (l + 4) = e 2S (l + ^). 



(4.10) 



(4.11) 



This junction condition gives the differential equation for the contour z(R) as seen from 
within the toroid which may be written 



4 = 4 + e ^- fe )-i (i+4). 



(4.12) 



For given V and rh both k(u) and z(u) are known on the toroid but z(u) and the constant k 



are as yet unknown. Explicitly k is given by (4.2) but with v given as a function of u and V 



by (4.8). It is useful in what follows to set 



1 + 4 , e> 



1 + 4. 



In terms of \ an d X the junction condition (4.11) may be written 



fe x = fe x or k + % = ^ + X- 



(4.13) 



(4.14) 



4.3 Evaluation of the "conicity" k 

k characterizes the conicity within the toroid. The conicity, defined as in [5], is the circum- 
ference of a circle of radius (R + dR) minus the circumference of a circle of radius R, that is 
2-Key dR divided by 2tt times the proper distance between the circles, that is 

= — = e k . (4.15) 

2Ke v + k dR 

Thus if A; > there is a deficit angle, a local "conicity" , typical of a circle on a cone centered 
on its apex or of spacetimes with a string. In the present case as we shall see there is 
"anti-conicity" , i.e. k < 0. 
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To calculate k we look at the point uo where z{R) reaches a maximum for a given V, 
that is where z 2 R = 0. Differentiating (4.12) with respect to u we find, writing a dash for 
derivatives with respect to u 



-znz'n = e 2 ^ [k'(l + z\) + z R z' R ] 



so where zr = 0, 



k' 



Z R Z R 
(1 + 4) 



(4.16) 



(4.17) 



Since k(u) and zr{u) are known, this equation may be used to determine the value of u, uo 



where zr = 0. We have here two conditions, one (4.17) that fixes uq and one from (4.14) or 



(4.12) that gives us k. Thus where u = uq, differentiating (4.14) where needed, 



x = x' = o 



k = k + x and k! + x' = 0. 



(4.18) 



From (4.8) we find that on the toroid 

v(u, V) = u 



2V 2 



and with v understood to be this function, we have, following (4.2) 



k = m 2 k(u,V) where k(u, V) = P(P\ -uP) + \v 



P 1 



(4.19) 



(4.20) 



Instead of seeking the values of k and uo where z(R) maximizes on a chosen potential surface 
V, we find it convenient to choose values of uq and V and to seek k and rh for which 



those values obey the maximizing condition (4.18). Using (4.20) we thus obtain rh and k as 



functions of uq and V from (4.18): 



k(uo, V) = m 2 k + x with m(uQ,V) 



-X' 



(4.21) 



We used Mathematica to make a parametric plot of k as a function of rh for different values 
of V. We see in Figure 3 that k is negative. For a given rh = 2irri , k becomes more negative 
with increasing V. 
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Figure 3: k as a function of in for, from top to bottom, V — 1.05, 1.10, 1.15, 1.20, 1.25, 1.35, 1.5. 

4.4 Surface pressures and surface mass-density on static toroidal shells 

We calculated the surface density and the surface stresses in the shell using the space-time 
fitting conditions of Israel. For a general position on the toroid these expressions prove 
tediously long and unenlightening when all the differentiations are inserted. They become 
somewhat simpler when v = e v = ±1 or sin 77 = i.e. on the equatorial plane where they 
reach extreme values. We obtain expressions for the energy per unit area and the tensions 
or pressures in the toroid: 



Notice that if one sets / = / = k = 1, the mass-energy density per unit length reduces to its 
classical value 



a, p v and p v are quite complicated parametrized functions of m and V. One possible test that 
the results are sensible consists in calculating the equilibrium of the forces in the equatorial 
plane in the limit of small rh. R is the radius of curvature in the cp direction in the classical 





kg = —2vdnip. 



(4.23) 
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limit. Let b be the curvature radius in the r/ direction. The non-relativistic equilibrium of 
forces at v = ±1 is given by 

P V Pr, _ 12 . W^llM) _ ! M n 4 x 

Now, 

1 fit — e v ) ,1 (u — e„) n— s r 

t = W and « = [i + 2( I ,-X)f] - e '^ T >- (425) 

We calculated with Mathematica the ratio in (4.24) for V = 1.2 and a mass m ~ 0.01. We 
found that for v = +1, the ratio is ~ 1.01 and for v = —1, it is ~ 1.02. These are reasonably 
close to 1. For a mass m ~ 0.2 which is rather relativistic the ratios are respectively ~ 1.2 
and ~ 1.6. cr± varies on both sides of the toroid in the equatorial plane in terms of the mass 




m 



Figure 4: The vertical dashed line is the limit where the energy condition (p v /a)— < 1 breaks down. 
Plain line: the matter densities or ana± as functions of m for V = 1.2. er_ is the lower curve. Thick 
dashed lines: (p v /c)± as functions of rh for V = 1.2. The lower line is for v = 1. The pressure 
is never too great. Tiny dotted lines: {p t p/a)± as functions of rh for V = 1.2. The lower line is for 
v = 1. For rh = 0.505 the pressure p v / 'a = 1. 

m. The characteristic property of <r_ is that it becomes negative for m > 0.75 beyond which 
the solution becomes unphysical. Figure 4 represents the ratio (p v /cr)± and can see that the 
pressure in the inside of the toroid, v = — 1 or r\ = it becomes too big for rh > 0.505. 
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1.0 1.1 1.2 1.3 1.4 1.5 1.6 



V 

Figure 5: The dominant energy condition limits the mass that can be put on the equipotential 
tp = V = X-. This maximum mass is plotted as a function of V . The limiting condition is (^r)- = 1. 
For the shape of the equipotentials see Figure 2. 

Figure 5 shows that for fh < 0.505 the pressure in the ry-direction is never too large. 
Thus for V = 1.2 all energy conditions hold for fh not too large: < m < 0.505. Figure 5 
represents the limits of fh as a function of V. The line is a polynomial interpolation. The 
toroids do not satisfy the energy condition above that line. 

Let a± , pr± and represent the energy densities and pressures in the equatorial plane 
at v = ±1 where we expect to find limiting conditions. The dominant energy conditions are 
a± > 0, — 1 < (pri/o')± < 1 and —1 < (p^/cr)± < 1. We shall illustrate the situation with 
numerical results which are typical of the general situation taking V = 1.2. 

5 A toroidal solenoid's metric and junction conditions 

In Section 3 we gave the general equatorially-symmetric static Weyl solution of Einstein's 
equations in toroidal coordinates. Our aim now is to fit this exterior solution to Bonnor's 
metric which has a toroidal gravomagnetic field. We fit on a torus of the exterior Weyl 
solution. This Weyl metric can be written with £ = e~^. At large distances r\ — > and 
u — > 1, \/u — v — > \piajr. Hence tp — > \/2(a/r) ]TJ ai = M/r. Notice that all the a/ contribute 
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to the asymptotic mass M, not just ao- We use the external potential (2.13) outside our 



toroidal shell where £ < Z that is u < u s = cosh Z. The same conicity problem occurs when 
we try to fit the primitive form of Bonnor's metric inside our torus. We therefore generalise 
Bonnor's metric to include an extra constant conicity term e 2k ; this is easily done since 
Bonnor's metric has a Killing vector and Einstein's differential equations are local. We may 
replace Tp in Bonnor's solution by any constant multiple of it and the metric will still be a 
solution locally. Of course the metric no longer satisfies the condition that it is regular on the 
axis. It will now have a string discontinuity there, but the axis is not included in the interior 
of our torus, so no such discontinuity occurs in the part of the space we use. Replacing Tp by 



e k ip in Bonnor's metric (1.8) and setting R = Re k we obtain the metric 



ds 



- - F[dt - (1 - F- l )dz} 2 - (e 2k dR 2 + R 2 d<p 2 + F~ x dz 2 ). 



(5.1) 



This is the one we use for the interior of our torus. Rewriting F as a function of R rather 



than R cf (1.9) ,F = 81 In (R/a) + C where C is another constant. To make easy comparison 



with our work on equipotential toroids in Section 4, we choose C = exp[— m^2u s P(u s )] so 
that the potential at R = a is that found earlier for the Bach-Weyl toroids. This formula 
defines our parameter m which is no longer the total mass which we call M. However for 
large V the Bach-Weyl toroids approximate tori so for them I = 0, m = M. Because of the 



change of metric from (1.8) to (5.1) we now have different field components, 'H (f = 8e I 



rather than 8/ cf (1.10). The constant k is determined by the same procedure as we used for 
static toroids. 



5.1 Fitting the potentials and the gamma metrics 



The internal metric (5.1 ) and the external metric (3.1 ) with the potential given by (2.13) and 



k by (3.12), must give the same induced metric on the torus. Comparing the coefficients of 



dt we find F = e~ 2 ^ hence from d2.13b we have 

1 



^2 aiPh{u s ) cos(Zry) = \n[F{R)]/y/u a - cos rj. 



Comparing coefficients of dip 2 we find 



R 2 = R 2 e 2 ^ = R 2 /F(R) = a 2 (u 2 - l)/[{u s - cosr]) 2 F}, 



(5.2) 



(5.3) 



which give R as an implicit function of rj on the torus. Thus F(R) in (5.1) is a 'known' 



function of 77, and the coefficients in the Fourier series of the right-hand side of (5.2) divided 
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by Pl{u s ) appropriately give us the a%. With those known the function k is known from the 
sum (3.12). Comparing coefficients of drj 2 we find, on re-ordering the terms and multiplying 
byF 



e 2k h 2 -Fe 2k R' 2 , 



(5.4) 



where a ' stands for a derivative with respect to 77. The unknowns are the function z(rj) that 
gives z as a function of r\ on the torus and the constant k; R(ri) is given as the solution 



to (5.3). The boundary condition is that z = at r\ = 0. We must still determine k; for 



any selected value of it we can imagine integrating (5.4) starting from -q = 0. Eventually the 



right-hand side will reach zero where z reaches its maximum. However in general d^ (right 



hand side) will not be zero when right-hand side = 0. But d ri (d v z) 2 = 2d„zd 2 z which must 



be zero when d„z = 0. Hence k must be so chosen that 



d{e 2k h 2 )/dr 1 = e" : d/dr ] [FR^} 



2k, 



(5.5) 



when 



e 2k h 2 



e 2k FR' 2 . 



Dividing (5.5) by (5.6) and letting rj t , pronounced eta-top, be the solution for rj of 



d v k = dr$Jx{F x l 2 l£ /h)] = d v ln 



R 



then evaluating (5.6) at £ = Z, r/ = rjt 



k=[k + ln(h/R) - i In F] z , vt ={k — ln[(R'/R) sinh^)]}^. 



(5.6) 



(5.7) 



(5.8) 



Thus we evaluate k and ensure that the two gamma metrics fit on the torus. As stated in the 



introduction equations (1.4), (1.5) and (1.6) constitute the complete set of Einstein's equa- 



tions for stationary space-times for which the Killing vector is time-like. As those equations 
hold both inside and outside matter and do not mention the vector potential A as opposed 
to the gravomagnetic field B, the boundary conditions implied by them for shell distributions 
of matter have that property too. The general fitting procedure of Israel involves the vector 
potential so it can be simplified for stationary metrics. These boundary conditions arise from 



integrating the Einstein equations (1.4), (1.5) and (1.6) across the bounding torus £ = Z. 
We integrate (1.4) and use the continuity of £ the coefficient of dt 2 to find the jump in the 



gradient of ip along the normal, denoting the integrated Too by tqo etc. Although there is a 
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step in the value of B 2 it does not itself have a delta-function so it does not contribute to the 
integral across the surface. 

- e~ 2 ^ J Vi;.dS = J Rov^gSx = k(t 00 - l - ? T )dS (5.9) 

so 

- e- 2 ^[n.W] = k[t 00 - g fr] = rea 00 . (5.10) 
This is the generalization of the electrical [n.E] = Air a. Evaluating the left-hand side 

e -fy-k h -ifop/Q£ + l e-^rixdhxF/dR = e^/T^F 3 / 2 sinhZip u - 4Je%/i?] = w o (5.11) 



where n\ = —z'/n and n = (e z') + -Fi? = (e h) . The integral form of (1.5) is 
J "Hdl = — 2re J JdS, where dl lies along the boundary of any chosen surface S. The only 
non-zero component is found by applying this to an elemental thin surface that cuts a small 
piece of the torus's surface orthogonally at constant rj, there is only a contribution to the line 
integral from inside because % is zero outside. So, remembering that x 3 = rj, 

fB^ = U 4> = 8Ie~ l = 2k^^R. (5.12) 

This is the exact analogue of the relationship between the discontinuity of the surface com- 
ponent of magnetic field and the surface current in electrodynamics. The line pressure in the 



(p direction is Kp v = ~f22T 22 which can be evaluated by integrating (1.6) through the surface 
of the torus. A similar calculation gives the sum 7fczT , but to evaluate both of these we 
need expressions for the integrals of the components of , the spatial curvature, that are 
transverse to the normal to the torus. Using Israel's method applied in the gamma 3-space 
these can be found from the external curvatures of the torus in the external and internal 
spaces between which it lies k%, k%. Following Israel for a, b transverse to the normals, (we 
take the normals to point into the volumes in which the external curvature is calculated, 
hence the result is a sum rather than a difference of re's). 



/ 



AVT^C = < + (5-13) 



In integrating (1.6) across the torus, the £ ;fc only has discontinuities along the normal since 
the potential is continuous, so £ ;fc; ' does not contribute to the purely transverse components 
of the delta function on the right-hand side. The contribution of the first term may be 
evaluated following Israel's method but one dimension lower since P kl is the Ricci tensor of 
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the 3-metric of 7-space. To apply Israel's formalism we calculate the external curvatures of 
the torus in the two gamma-spaces, with the normals pointing into each space in turn, and 
then add the results. We take 9 a = (ip, rj), a = 2, 3 to be the coordinates on the torus itself. 
Outside: 

The normal in the external space is along — V£ so n k = (— e^ +k h, 0, 0) in the gamma- 
metric '■)kidx k dx l = e 2 ^(e 2k h 2 d( 2 + R 2 d(f) 2 + e 2k h 2 dn 2 ). The external curvature is 

dx k dx l dx k dx k dx l 

Kab = ~WW nk ' 1 = ~W bTlk + Wd¥> klTlm ' ( ' 

where the are the affine connections of the gamma-space. Since the normal has only a 
first component the first of the two terms on the right-hand side vanishes and we need only 
calculate the \\ t with k, I = 2, 3. The only surviving terms of this type are, 

\\ 2 = -\ e- 2 ^- 2k hr 2 d^R 2 ) ■ Xl 3 = -d c (i; + k + lnh)/dC. (5.15) 



From these we deduce k?, = X\ 2 n\e 2 ^/R 2 and n 3 = Aj> 3 nie 2 '^ 2k /ti 



2 



4 = e-^v^Tf^K ~v) + fa^]; (5-16) 
4 = e-^- k ^u 2 s -l[{^ u + k u ){u s -v)- 1]. (5.17) 

Inside: 

Within the torus the gamma-metric is j kl dx k dx l = e 2k dR 2 + R 2 dcp 2 + F~ 1 dz 2 , where 
x k = (R, if, z) and on the torus it is R 2 dcp 2 + Hdrj 2 , where H = e 2k R 12 + z' 2 /F. The normal 
that points into the torus is n k = (ni, 0, 7x3) = (— z', 0, R')/n. Where n is the quantity 
defined under (??) and 

_ dx k _ _ dx k dx l ^m_ . , n , 

Kab = ~W bUk + d¥d¥> klUm - ( ] 

i 3 

As n k has components in both the 1 and 3 directions, we now need both \ kl and X ki with 
k, I = 2, 3. However 7 n is now constant and both 7 22 and 733 depend only on R therefore 
= \\ t = A23 = A33 = A~2 2 = X 3 n . We need 

4 = -e~ 2k R; A 33 = - e" 2 ^^; A? 3 = A = ~~ ^ (5-19) 
22 33 2 F 2 dR 13 31 2 dR V ' 

Thus we find K22 = A^ni; ^33 = — (R'n'i + z'n' 3 ) + 2R'z' \i 3 n 3 + z /2 A 33 ni. 

Now —(R'n^ + z'n 3 ) = (z' /R')'(R') 2 /n so evaluating k 2 and k 3 we obtain, 



Rn 



^) {Ry + 2R"z'X[ 3 -z! i X m 
R / 
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jn. (5.20) 



The integral across the surface of the sum of the transverse spatial curvatures is found from 



(1.6) 



«i + K\ + 4 + «i = K{ lkl T kl + T) = ^. (5.21) 



So tqo is determined and hence r is known from (5.9). Using the equation above, the integral 



of the 22 component of (1.6) yields on multiplication by 722 

The surface energy density, a and the p^ principal components of the surface stress, are 
related to the components of the surface energy tensor through the relationship 

T = a — ptp — p v (5.23) 

and those that depend on the velocity v which is in the 77 direction, 

too _ <r + p v _ n . t[_ (a+p v )v , , 

With p,p known, these three may be solved for v ,p v , and a successively giving, 

W t3/£ 

■ W= - ) ■ -, (5.25) 



"1 + ^13^72 ' roo/e 2 - l(r + Pip ) 

and 

(l-^ 2 )(Too/e 2 )-(r+p y ) 
Pr?= l^v^ o- = T + ^+p^. (5.26) 

The dominant energy conditions are a > \p v \ and <r > \p v \. 

6 Exploration of the Rolling Tori 

6.1 A relativistic example 

Apart from an overall scaling there is a three dimensional set of solutions. We take the radius 
of the line torus as our unit of length, GM / c 2 gives another length while the total mass flux 
/ is a mass per unit time which becomes dimensionless when we use ct for the time. We 
find it convenient to use not the total mass M itself but our quantity m that is more closely 
related to the potential inside the torus. In practice M is within 10% of m. There will be 
limits on the size of both the mass and the mass flux but these we shall explore. Since the 
radius of the line torus is our unit of length, the mass will be limited by the condition that it 
is supported by pressure and does not collapse into a black hole. Likewise a large mass flux 
can only be held in by a torus of small minor diameter. Given the two parameters m and / 
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we may still choose the minor axis on the torus on whose surface the mass moves. It can be 
small giving us a narrow tube or large like a bulky cored apple. Again there are limits on the 
minor axis caused by the energy conditions on the principal stresses and the surface density 
which must be positive. These imply that the velocity v/c does not exceed unity. 



J7/JT 




0.0 0.2 0.4 0.6 0.8 1.0 



Figure 6: Properties of the torus with m/a = 0.659, / = 0.008, b/a — 0.297. The rolling velocity rises 
to 0.85c as i] increases from to tt, but the positive surface density decreases over that range. This 
system is strongly relativistic. 

Figures 6 and 7 describe a strongly relativistic torus with m/a = 0.70, I = 0.008 and 
a minor semi-axis of (I — u 2 s )~ l / 2 a = 0.297. The £ on this torus is 1.93. Figure 6 gives the 
velocity v in the rj direction as a function of rj. It gets to 0.85 at the inner equator but it 
must be emphasized that this velocity is not along the equator but at right angles to it in the 
T) direction. This maximum occurs on the inner equator but this coincides with the minimum 
surface density, a, that is also portrayed in figure 6. It is positive as it should be. Figure 7(a) 
gives the ratios p^/cr, p„/ a which must lie between one and minus one if the dominant energy 
condition is to be satisfied. p v /a reaches 0.934 at the inner equator so it is only just below 
the limit showing that we have a highly relativistic system. Notice that p„ starts positive but 
declines through zero at r\ = 3ir/4 whereafter it becomes significantly negative showing that 
tension is necessary to supplement gravity in opposing the centrifugal force due to the fast 
rolling motion near the inner equator. We now give the other parameters of this particular 
system. Outer and inner equatorial radii a a /a = 1.21; ai/a = 0.826. Fourier coefficients 
a* = {0.445, 0.0205, 4.44 x 10" 4 , 9.32 x 10" 6 }, % = 0, 3, the true total mass M/a = 0.659 
as compared with m/a = 0.7. Notice that this is a higher total mass than the maximum that 
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Figure 7: The energy conditions for the same torus as Figure 6. The ratios of the principal pressures 
to the surface density as functions of 77 must remain between ±1. These are satisfied but only just as 
the maximum is 0.934. Notice that p v becomes negative, (a tension) helping gravity to balance the 
high centrifugal force where the rolling velocity is large. 

can be supported as a Bach-Weyl toroid, or as a static torus (see figure 8). The conicity is 
k = +0.0192 and z achieves its maximum at 77 = 0.361-7T. 

6.2 Tests of the computation in the classical limits 

Of course our solutions which are worked out via the relativistic computation of k also contain 
classical cases when the potential is everywhere small and the velocities much less than 1. 
These provide useful checks on our numerics because in the static case the pressures times 
the external curvatures in their directions must balance the gravitational pull on the density. 
This balance gives 

test : (kIp^/ct + re|p r? /cr)/(recr/4) = 1 (6.1) 

This test is satisfied to an accuracy of better than half a percent for all rj on a static torus 
with rh = 0.001. A similar test but with a velocity is also satisfied when the static test is 
strongly violated for this classical system whose velocities exceed y/p/a. Here we must use 
our dynamic vtest which is, 

vtest : {rfiptp/a + K\{jp n /a + i> 2 ))/(ko-/4) = 1 (6.2) 

This test on a system with fh = 0.001 and / = 0.0001 with a maximum velocity of v/c = 
0.025 is satisfied to a similar accuracy even when the static test (applied wrongly to this 
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dynamic case) is strongly violated. Evidently inclusion of the centrifugal force makes a very 
considerable difference. 




Figure 8: Contours at m/a = 0.25,0.5,0.6,0.7,0.75 of the limiting surface in the I,b/a,m/a space. 
Within the contours there are solutions obeying the dominant energy condition. Notice that the higher 
contours overhang the lower ones at the upper right where the lower ones are dashed. In most places 
the limit comes from p v /<7 = 1 except at the upper left where it comes from p v /cr = —1. Symbols 
denote computed points. 

6.3 Relativistic limits due to the energy conditions 

Figure 8 shows the region in b/a, I, m/a-space within which the solutions satisfy the dominant 
energy conditions. The contours are at constant m/a. Evidently the rolling motion allows 
greater values of m/a than those available for static toroids. Notice that the contours at 
larger m/a overhang at the upper right those at lower m/a so the contours actually cross. 
We have not seen a topographical map where this happens but wherever a mountain has 
been undercut by a glacier leaving an overhang it should be expected. Where the lower 
contours are overlaid by higher ones they are dotted. Over most of the diagr am it is the 
Pip/p = 1 that determines the limiting surface but this is replaced by p^/o = —1 at small 
b/a and large I. In all cases the surface density falls as r/ increases so the highest surface 
densities are achieved on the outside equator. Generally the velocity increases inwards but 
for thin tori this can go the other way and there are even cases in which the maximum in the 
velocity is not on the equatorial plane. When the velocities are not very high both principal 
pressures are p sttositive but at high rolling velocities p v becomes negative, so tension is 
required as well as gravity to hold the matter in against the centrifugal force of the rolling 
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motion. We decided not to burden the reader with a full panoply of different cases but to 
give figure 8 whose determination required many solutions to be run in the neighbourhood of 
this bounding surface. At the bottom of figure 8 the static tori at fixed m/a do not extend to 
b/a = 0. This again demonstrates that static line tori disobey the energy conditions. Looking 
a bit higher up we see this is also the case for our rolling tori, in both cases p^/cr becomes too 
large. In choosing Bonnor's metric inside our torus we imposed a structure on the internal 
potential that is not that which would be found if the torus were cut and unrolled into a 
cylinder. Thus there will be many more solutions for rolling tori that are not included in the 
particular set we have chose stn. However the simplicity of Bonnor's solution suggests that 
those investigated here will be among the simplest solutions of this type. 

Acknowledgments 

We thank Jin Bicak for his interest and advice and John Harper for a reference. We par- 
ticularly admired the meticulous referee who rightly rejected our earlier paper on the statics 
only. Nevertheless the Mathematica computations there (reproduced here) were correct. 
We aspire to meet his high standard of accuracy and typography in this more complicated 
problem. 

References 

[1] Olver F W J, Lozier D W, Boisvert R F and Clark C W 2010 NIST Handbook of 
Mathematical functions formula 14.8.14 (Cambridge: CUP) 

[2] Bicak J, Ledvinka T, Schmidt B G and Zofka M 2004 Static fluid cylinders and their 
fields: global solutions Class. Quantum Grav. 21 1583 ( Preprint |arXiv:gr-qc/ 0403012) 

[3] Bach R Weyl H 1922 Neue Losungen der Einsteinschen Gravitationsgleichungen Math. 
Z. 13 134 

[4] Bonnor W B 1969 The gravitational field of light Comm. Math. Phys. 13 163 

[5] Frolov V P, Israel W and Unruh W G 1989 The gravitational field of straight and 
circular cosmic strings: Relation between gravitational mass, angular deficit, and internal 
structure Phys. Rev. D 39 1084 



26 



[6] Geroch R 1971 A method for generating solutions of Einstein's equations J. Math. Phys. 
12 918 

[7] Hoenselaers C 1995 The Weyl solution for a ring in a homogeneous field Class. Quantum 
Grav. 12 141 

[8] Landau L and Lifshitz E 1966 Theorie du champ (Moscow: Mir) pp 354, p355 §95 
Problem 

[9] Morse P M and Feshbach H 1953 Methods of theoretical physics (New York: McGraw- 
Hill) 

[10] Semerak O, Zellerin T and Zacek M 1999 The structure of superposed Weyl fields Mon. 
Not. Roy. Ast. Soc. 308 69 



Appendix 

A Evaluation of the indefinite integrals involved in equa- 



tion (3.11) 



Our general solution for k was given in the unwieldy form (3.11 ) in terms of indefinite integrals 



of polynomials times products of Legendre functions. Here we show how it may be expressed 
in terms of Legendre functions themselves. The recurrence relation (u 2 - 1)P' L = L(uP L - 
Pi_ i) may be used to eliminate derivatives with respect to u in the final results. 

We first show that all eight integrals may be reduced to known functions together with 
j PiP^jdu and j PLPMudu and then evaluate those integrals. The second, fifth, seventh and 



eighth terms on the right of (3.11 ) are already of the above form. The third is, integrated by 
parts, 

/ (PlPm)'(u 2 - l)du = P L P M (u 2 - 1) - 2 j P L P M udu (A.l) 

which is of the desired form. The first term gives two integrals both of which are l,m 
symmetric. Integrating by parts and using ( |A.l ) and Legendre's equation ( |2.11 ), 



J P' L P' M (u 2 ~ l)udu = P L P' M (u 2 - l)u - (m 2 - § J P L P M udu - J P L P' M ( 



l)du 



\(PlPm)'(u 2 - l)u - \ (I 2 + m 2 - \ ) J P L P M udu - 1 J (P L P M )'(u 2 - l)du 
\ {u 2 - 1) [(PlPm)'u - P L P M ] + [1 - \ {I 2 + m 2 - \ )] J P L P M udu, 

(A.2) 
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which is of the desired form, and 

J P' L P' M (u 2 - l)du = P L P' M (u 2 - 1) - (m 2 - \) f P L P M du 

= | (PlPm)'(u 2 - 1) - \ (I 2 + m 2 - \ ) J P L P M du. (A.3) 

The fourth and sixth term get the required form as follows: 

J P' L P M (u 2 - l)du = P L P M (u 2 - 1) - 2 J P L P M udu - j P L P' M (u 2 ~ l)du, (A.4) 
which may be written 

2 J P' L P M (u 2 -l)du = P L P M (u 2 -l)-2 J P L P M udu- J (P L P' M -P' L P M )(u 2 -l)du (A.5) 
The last term can be simplified because from Legendre's equation for Pl and Pm, 

[(u 2 - 1)(P L P' M ~ P'lPm)}' = -{I 2 ~ m 2 )P L P M . (A.6) 

Thus, 

J (PlP'm ~ P'lPm)(u 2 - l)du = (P L P' M - P' L P M )u(u 2 - 1) + {I 2 - m 2 ) J P L P M udu. (A.7) 

2 J P' L P M {u 2 - l)du = (u 2 - 1)[PlPm- (P L P' M - P' L P M ) u] - J P L PMudu. 

S A A S 

(A.8) 

Where symmetrical and anti-symmetrical terms for l,m interchange are indicated. Taking 
account of symmetry and anti-symmetry in the coefficients of m in the 4th and 6th terms 



of terms of (3.11) we symmetrize for the l,m interchange and obtain their contribution to 



the coefficients of a;a m . To evaluate f PLPMdu we integrate Pm times equation (2.11) and 
subtract the result with L and M interchanged to obtain 

(l 2 -m 2 ) J P L P M du = -( U 2 -1)(P L P' M -P' L P M ) = {l-m)uP L P M -LP L ^P M + MP L P M -i. 

(A.9) 

This gives us the desired integral whenever I ^ m. The I = m integral is given by differenti- 
ating with respect to L and putting M = L: 

21 J P 2 du = uPl - P L ^P L -(I -I) (PL d -^ ~ ^Pl-i) ■ (A.10) 
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To evaluate J PlPmucIu we add (u 2 — 1)P' l Pm and (u 2 — 1)P' m Pl , integrate by parts, and 
use the recurrence relations to obtain 

(l + m+1) J P L P M udu = P L P M (u 2 - 1) + (I - i ) J P L -iP M du + (m - § ) / P L P M -idu. 

(A.ll) 

Hence, 



PlPmucIu 



1 

l + m + 1 



(u 2 - 1)P L P M + j Ph-xPudu + {m-\) J P L P M -idu 

(A.12 



which are evaluated in (A. 9) and (A. 10). 



The different terms in equation (3.11) have different dependencies on r\. We therefore 



collect the terms in cos (nrj) for different values of n involved. 

Collecting our results k can be written in terms of kf m (u) with n = ±1 or 0: 



k 



1 E£o Em=o a i a ™ c(l + m + l)kj + c(l + m)kf +c(l + m- l)fc, i 



+ 



i 1 v—voo v— VOO 

+ 8 2^i=o 2^ 



m =0 a l a m 



c(l-m+ l)k}_ m + c(l - m)kl_ m + c(l-m- 1)^1 (A.13) 



where 



1) [(PlPm)'u + (l + m)P L P M + {l- m){P L P' M - P' L P M )u] 



udu. 



| [(/ _ ra)a _i]( i + ra _i)} J P L P 1 

/U 
P L P M du 



h u — 9 



1) [(PlPm)'u -(1 + m)P L P M -{I- m){P L P' M - P' L P M ) 



u 



+ 



{[{I 



m 



l](l + m + l)\ / P L P M udu, 



(A.14) 
(A.15) 

(A.16) 



The coefficients below have -m written for m in the above in leaving L and M unchanged; 

fcj-m = ~(u 2 - 1) [(PlPm)'u + {l- m)P L P M + (I + m){P L P' M - P' L P M )u] 

- {[(/ + m) 2 -l](Z-m-l)}^ P L P M udu. (A.17) 

fc£_ m = 2[(P L P M )'(u 2 - 1) - (I + m) 2 J P L P M du]. (A.18) 
Notice that when I = m we need the final integral except when I or m is zero. 



v l,-m 



-{u 2 - 1) [(PlPm)'u -(I- m)P L P M - (I + m)(P L P' M - P^P, 



M)U 



{[(1 



+ < (l + m) 



\}{l-m + l) \ J P L P M udu, 



(A.19) 



29 



which completes our calculation of k? which have been checked numerically via numerical 



integration of equation (3.9) at constant r\ and Fourier transformation. For the rolling tori 
the aj fall rapidly as i increases so only a few terms are needed for accuracies of 1 part in 
10,000. 
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